Field-induced domain wall propagation: beyond the one-dimensional model 
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We have investigated numerically the field-driven propagation of perpendicularly magnetized fer- 
romagnetic layers. It was then compared to the historical one-dimensional domain wall (DW) 
propagation model widely used in spintronics studies of magnetic nanostructures. In the particular 
regime of layer thickness (h) of the order of the exchange length, anomalous velocity peaks appear 
in the precessional regime, their shape and position shifting with h. This has also been observed 
experimentally. Analyses of the simulations show a distinct correlation between the curvature of the 
DW and the twist of the magnetization vector within it, and the velocity peak. Associating a phe- 
nomenological description of this twist with a four-coordinate DW propagation model, we reproduce 
very well these kinks and show that they result from the torque exerted by the stray field created 
by the domains on the twisted magnetization. The position of the peaks is well predicted from the 
DW's first flexural mode frequency, and depends strongly on the layer thickness. Comparison of the 
proposed model to DW propagation data obtained on dilute semiconductor ferromagnets GaMnAs 
and GaMnAsP sheds light on the origin of the measured peaks. 



PACS numbers: 75.50.Pp,75.60.Ch,75.78.Fg 
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I. INTRODUCTION 

The 1970's saw the establishment of the main the- 
ories describing magnetic domain wall (DW) propaga- 
tion, among which the so-called one-dimensional (ID) 
model 1 - - — . These studies fueled the intense efforts towards 
the building of magnetic bubble memories. They have 
recently come back into fashion with new schemes being 
proposed to use domains or domain walls as the building 
blocks of a non-volatile and down-scalable memory^—. 
Described in different formalisms 1 - - — , the ID model as- 
sumes that the DW propagation can be fully described 
by two time-dependent variables: the position go(^) and 
azimuthal angle 4>o{t) of the magnetization unit vector m. 
For a defect-free sample, the velocity versus field curve 
v(H) is then shown to consist of a high mobility regime 
where the configuration of rh within the DW is stationary 
(</>o(i) constant) up to the Walker field Hw, followed by 
a linear, lower mobility regime where fh precesses around 
the applied field. 

While these two regimes have indeed been evidenced 
experimentally^—, one or two unexpected kinks have re- 
peatedly been observed in a variety of configurations: 
in-plane magnetized Permalloy^, out-of-plane magne- 
tized garnets^ 3 , as well as out-of-plane magnetized fer- 
romagnetic semiconductor GaMnAs 11 ^. It has also been 
observed in field-assisted current-induced propagation 14 . 
The present work aims to address this issue, all the more 
crucial as the ID model is now routinely used as a guide- 
line for a vast number of studies on current-induced DW 
propagation 1 ^—. Whereas differences with the ID model 
in the stationary regime were attributed quite early on to 
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FIG. 1. [color on-line] (a) Geometry of the simulation, side 
view of the l(jmx/i layer. The DW propagates towards the 
left, its mean position along y is given by go (t) (not to scale), 
(b) DW propagation evidencing curvature, shown here for 
h=80 nm to highlight the effect, (c) DW free oscillation 
(H z =0), after launching from either of the top two positions. 



the nucleation/annihilation of Bloch lines^, those in the 
precessional regimes have, to our knowledge, only been 
considered by the numerical simulations of Patterson et. 
at^ who did not isolate the physical mechanism(s) for 
these kinks. 

In the following section of this paper (Section II), 
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we present simulated v(H) curves for perpendicularly 
magnetized ferromagnetic layers of increasing thickness 
(/i=10-40 nm) obtained with an open source micro- 
magnetic simulation software 2 ^. In Sections III and IV, 
we analyze extensively the three component magnetic 
configuration within the DWs to relate it to the position 
of the v(H) peaks. In particular, we evidence the role of 
the curvature of the DW and the twist of its magnetiza- 
tion vector in increasing the DW velocity above the ID 
model value for certain fields. Those particular fields are 
found to correspond to the first flexural mode resonance 
of the domain wall. Finally, in section V this approach 
is tested against experimental v(H) curves obtained on 
GaMnAs and GaMnAsP thin films of varying thickness, 
magnetization, magnetic anisotropy and exchange con- 
stants. Ferromagnetic semiconductors are ideal for these 
studies as their weak pinning allows all DW propaga- 
tion regimes to be observed^, their magnetic parame- 
ters can easily be tuned by careful dopingSir— , and the 
layer thickness can be varied across the exchange length. 
Section VI summarizes the results and opens new per- 
spectives for DW propagation in other geometries. 



II. SIMULATIONS 



A. Geometry &: results 

The simulation was laid out on a 2D mesh with cell 
dimensions c 2 , and a strip of length 2 - d=l fim. Vary- 
ing heights were investigated: h=\0 nm (c=l nm), and 
h—20, 30, 40 nm (c=2 nm). The film was taken in- 
finite in the x direction (Fig. QJi). The applied field 
H=H Z was constant and perpendicular to the plane 
of the layer, and micromagnetic parameters typical of 
GaMnAs 2 ^ were used: magnetization M=33 kA m -1 , 
uniaxial anisotropy coefficient if M =8878 J m -3 and field 
HqH u = 2K U /M, exchange constant yl=5.10~ 14 J m _1 , 
and damping coefficient a=0.3. Given the high value of 
the factor Q — -^jjjfe >3 in our layers, the complex in- 
plane anisotropy of GaMnAs was ignored. Its influence 
will be discussed later. The magnetization orientation 
vector rh was defined in Cartesian (m x , m yi m z ) or (0, ip) 
coordinates (Fig. QJ,) . The position coordinate qo is given 
as usual 5 by 0=-^- sin where the dot denotes the time 
derivative. The instantaneous velocity Vi ns t = f 9< m' > 
was averaged over many periods and the resulting v(H) 
curves are shown in Fig. [2 For comparison, the baseline 
of the curves was fitted to the ID model velocity (keep- 
ing a fixed to 0.3, and varying {jlqHw and A), using the 
(go,0o)(t) solutions to the Landau-Lifshitz-Gilbert (LLG) 
equation 1 for H > Hw- 
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In these expressions, Hy/ = aM/2 is the Walker field 
for a film of infinite thickness, A=^ A/K the static DW 



width and T 



the precession period 



M07 y/m-Hl, 

(time for c6 (i) to span 360°). 

The simulations were done up to 300 mT, but only 
the first part of the curve is shown in Fig. [2] as the ve- 
locity maintains the same mobility for higher fields. For 
h=10 nm (Fig[2^), the numerical simulations as well as 
the ID model lead to a very standard curve. A high mo- 
bility is observed up to the Walker field, then a narrow 
negative mobility region followed by a positive lower mo- 
bility regime. The DW velocity reaches over 30 m s" 1 at 
/zo-ff=250 mT. No kinks were observed in the 0-300 mT 
range. The h— 20-40 nm curves evidence a similar behav- 
ior up to Walker breakdown, but distinctive kinks appear 
above: a broad peak («70 mT) for h=20 nm, a slightly 
asymetric peak (s=s37 mT) for h=30 nm, and two close 
peaks (25 and 35 mT) for h=A0 nm. Walker breakdown 
occurs at low field (fioHw=3 mT, velocity Vw=4 ms ) 
for the thinnest layer, and then gradually converges 
to the infinite thickness film value /j,qHw=G mT for 
h=40 nm. In the mean-time, A decreases slightly with 
thickness (from 2.6 nm to 2.3 nm), reflecting the change 
of demagnetizing factors of the DW with layer thickness, 
which in turn modifies the values of Hw and A3§. The 
precession frequency of the magnetization f pre c obtained 
as the inverse of the period of m x (t) was found identical 
for all thicknesses and linear with the applied field (solid 
line in Fig. [5]). 



B. Image analysis 

The three component m vector field images (of the type 
in Fig. [T}d,c) were then treated numerically to extract for 
each iteration the domain wall width A and the depth- 
dependent DW position q(z,t) and ip(z,t) angle (aver- 
aged over 2 A along y). The z-averaged angle <j>(t) is given 
by i ^2i' i>i[z, t) with Ni the number of cells across the 
thickness. Plotting simultaneously Vi ns t(t) and (j)(t) for 
H > H\y over one precession period shows the main fea- 
tures of the precessional regime (Fig. [3^,): m precesses 
around the applied field inducing a demagnetizing field 
across the wall. This sinusoidal field produces a torque 
on m leading to an alternating forwards or backwards 
motion of the DW. The applied field term in Eq. (1) 
insures that the total wall displacement over one period 
is strictly positive. As expected from the ID expres- 
sion for <po(t) (Eqs. (1,2)), the instantaneous velocity is 
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FIG. 2. [color on-line] Simulated velocity versus field 
(squares): (a) h—10 nm, (b) h—20 nm, (c) h—30 nm, and 
(d) h=40 nm. In (a), (b) and (d), the solid line is a fit to 
the ID model, and the stars the maximum DW elongation 
r)o. In (c) the open circles are the twist amplitude So and 
its Lorentzian fit (solid line). Investigated thicknesses he on 
either side of the exchange length, nA—27 nm. 



then maximum for </>=45° modulo 180° and minimum 
for </>=135° [180°], corresponding to the extrema of the 
demagnetizing torque (arrows in Fig. . 

Meanwhile, the DW length is not constant: the DW 
undergoes a T-periodic flexion (Fig. [Tb and|3b). Let us 
define the length of the domain wall l(t) as the curvilinear 
integral of q(z, t) and its elongation as 77(f) = [/(f) — h]/h. 
The simulations show that the DW goes from no elonga- 
tion (l=h, 77=0) to the maximum 770 while <p(t) varies 
from to 360°. Careful analysis shows that the curva- 
ture is always of an 77=1 mode type (following the num- 
bering of Slonczewski22), i.e odd about the layer's mid- 
height. For the thinnest h—10 nm layer, the amplitude of 
elongation increases progressively with field up to 1.53% 
for 300 mT (Fig. Hi). For ft,=20-40 nm, 770(H) reaches 
a maximum at a resonance field H res which decreases 
with layer thickness. This maximum elongation (stars 
in Fig. [2}j,c,d) increases with h: rjo—2.3% for h=20 nm 
(fiQH res —70 mT), 3.7% for /i=30 nm (/jioH res =35 mT) 
and 4.8% for h=i0 nm ([j, H res =25 mT). Far from H res , 
the flexion becomes negligible. Note that the H res fields 
correspond very exactly to those of the kinks in the v(H) 
curve. The enhancement of velocity therefore seems re- 
lated to the field-dependent elongation of the DW. 

For ft=10-30 nm and fi=40 nm when /iqH <30 mT, 
the twist amplitude is minimum for f=T/4, maximum 
around 3T/4, and null for f=0 [T/2] (Fig. [5b). Twist 
and curvature are intimately related, being in quadra- 
ture phase with respect to each other. The depth depen- 
dence ip(z,t) of the propagating DW evidences a clear 
twist from top to bottom of the layer (Fig. 0J, but quite 
differently from the text-book 180° twist case of a static 
Neel-Bloch-Neel Plotting ip(z,t) over a whole 

period shows that the twist is sinusoidal across the layer's 
thickness, and that its amplitude S(t) also varies sinu- 
soidally in time over T with a span of So (Fig. 2^). The 
time and depth dependence of tp(z,t) are then described 
in first approximation by the following phenomenological 
expressions: 



ij)(z,t) = (t> {t) + i)(z,t) 

~ T! 

ip(z,t) = -S(t,H) cos 



S(t,H) = -S (H)sm 



h 
2irt 

T{H) 



(3) 
(4) 

(5) 



The thickness dependence of Sq is as expected: it in- 
creases with thickness, since the exchange energy cost of a 
twist can be more easily accommodated in a thicker layer. 
Strikingly, Sq(H) follows the same trend as t]q(H): it 
reaches a maximum at H res (see Fig. [2]: for the fi=30 nm 
case). At resonance, the full twist amplitude of the DW 
magnetization across the layer's thickness (2So) reaches 
about 90° for h=20 nm, 216° for h=30 nm and 230° for 
h=40 nm. 

For the thickest layer /i=40 nm, a different behavior 
appears around [ioH=30 mT, as can be seen in Fig. \^p: 
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time 



the twisting direction inverts at a given depth. This cor- 
responds to a horizontal Bloch line (HBL). These appear 
where the DW is most labile as a result of the competition 
between the stray field emanating from the domains, and 
the demagnetizing field within the DW— . For h~30 nm, 
very loose HBLs also appear around 45 mT. This is con- 
sistent with the estimated HBL width ttA=tt^/Q~27 nm. 
A 40 nm layer can thus accommodate with difficulty a full 
HBL. 

To summarize, this first analysis gives two fundamen- 
tal results: the kinks in the v(H) curves occur at fields 
maximizing both the curvature and the twist of the DW, 
and twist and curvature are coupled variables. 



III. ANALYSIS 

In order to get more insight into the respective roles 
of the DW flexion and magnetization twist on the kinks 
of the v(H) curves, the analytical model of Malozemoff 
and Slonczewski2£ will be used. Let us beforehand re- 
mind the general equation coupling the time derivatives 
of (<Z0j<A))(t) in the ID model: 



FIG. 3. [color on-line] h=20 nm, fj,oH=70 mT (a) Instanta- 
neous velocity (open squares) and mean </>(t) (full squares) 
over one precession period. Arrows indicate maximum (mini- 
mum) velocities for the canting angle close to 45° [180°] (135° 
[180°]). (b) Domain wall elongation (solid line) and domain 
wall twist S. 
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FIG. 4. [color on-line] Depth dependence of %j) (z,t) over half 
a precession period for (a) ft=20 nm, /j,oH=70 mT and (b) 
h=40 nm, /io-ff=35 mT. Loose horizontal Bloch lines appear 
in the thicker layer. 



4o(t) 



7A sin 20o + aA(j) (t) 
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In the precessional regime, the propagation is achieved 
by a balance of the torque associated to the oscillating 
demagnetizing field (first term of the right hand side of 
Eq. (6)) created by the precession, and the dissipation of 
the energy brought by the field in the form of a damping 
contribution (second term). For H » Hw, the demag- 
netization term averages out to zero, and the damping 
term takes over and propels the DW forward with a ve- 
locity proportional to the applied field. In this simple 
ID model, two different routes appear to increase v(H): 
an increased precession rate or an increased damping. 
As mentioned above however, the simulations showed no 
substantial difference for 4> between the different thick- 
nesses. Concerning the damping, an alternative way for 
the system to dissipate energy would be by direct cou- 
pling to spin-waves (SW) in the bulk of the material, 
or by a spin-wave like excitation of the magnetization 



in the vicinity of the DW 



Plotting the 



X :V: Z 



(y) 



components for different thicknesses does not evidence 
any particular wake, but only a slight deformation of the 
DW profile for H w H res . Concerning bulk spin- waves, 
their lowest (fc=0) frequencies^ lie much higher than the 
typical precession frequencies, thus forbidding any direct 
coupling between the DW and bulk SWs. It therefore 
seems necessary to go beyond the ID model, and take 
into account the full thickness dependence of the magne- 
tization configuration. 

For a twisted domain wall taken infinite in the x direc- 
tion, the DW surface energy is expressed to first order in 
exchange, magnetostatic, and curvature terms as^: 



5 



2AA (W) 



(7) 



fi M 2 Asm 2 ip - Afj, TT MH y sin ip - 2/j, MHq 



ID DW sur- 
corresponds to the 



In this expression, <tq = 4 V AK is the 

face energy. The term 1 + \ (Vq) 2 

energy increase due to the elongation of the DW, and 
the term A (V^) 2 corresponds to the exchange energy 
cost arising from the twist. By analogy with a spring 
or a twisted rope, and looking at Fig. [3b, the DW en- 
ergy contains a kinetic (flexion V 2 q) and an elastic (twist 
V 2 t/>) term which balance each other out. When the DW 
is very twisted, the energy cost of an additional curvature 
is too high: the flexion is minimal, and vice- versa. H y can 
stem from volume or surface magnetic charges created by 
the DW or domains. Here, we will only consider H y as 
the y component of the stray field arising from up- and 
down-domains on either side of the DW. The curvature 
of the domain will be considered small enough to neglect 
the fields along y and z resulting from magnetic charges 
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appearing along the DW due to its curvature 2 ^. Among 
the different expressions demonstrated for Hy£>2&, we will 
use that of Ref. 0, as it remains valid down to the A « h 
limit. 



z 2 + A 2 /4 
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z arctan 



(h-z) 2 + A 2 /4 
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It is determined by computing the potential arising 
from the -M and +M surface charges (varying linearly 
over a width A) of domains on either side of the DW. Any 
effect of the volume charges created within the DW are 
ignored. This field is stationary, but varies throughout 
the thickness of the material. It pins the magnetization 
of a static DW in Neel configuration at the surfaces of 
the layer, so that the equilibrium structure of the wall 
is twisted. For h—30 nm, it amounts to a field of over 
55 mT at the surface of the layer. This is likely an over- 
estimation of the actual H y though, since the volume 
contribution has been ignored. Finally, using the LLG 
equation leads to the equation^: 



q(z,t) 



7A 
1 + a 2 



MoM 7T 2AV 2 V> aa W 2 g 
— ^— sm2</> + an H - -n Q H y cosip ^— + — 



(9) 



The DW velocity is then obtained by a double integra- 
tion on time and space v(H)=^ J t _ Q ^ J z=0 q(z,t)dz dt, 
provided the canting angle is known. As expressed in 
Eqs. (3)-(5), tp{z,t) can easily be written as the sum 
of the ID model uniform </>o(t) given by Eq. ([2]), and a 
depth-dependent angle ip(z,t) varying sinusoidally across 
the layer (in the absence of HBL). Moreover, the oscil- 
lation amplitude Sq{H) can be fitted by a Lorentzian 
(Fig. [2J:, solid line). The negative signs in Eqs. (4,5) 
mimic the phase seen in the simulations: the twist is null 
when the wall is in a Bloch configuration (cj)o—0 [180°]), 
and extremum when it is in a pseudo Neel configuration 
00=90° [180°]). 



The resulting computed velocity is plotted in Fig. [5]for 
/i=30 nm, along with the corresponding micromagnetic 
simulation already presented in Fig. The Walker field 
and DW thickness injected in Eq. (2) are those of the ID 
model baseline fit of the h—30 nm curve : fj,oHw=5 m T 
and A=2.28 nm. The v(H) curve obtained using the 
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analytical model given by Eqs. (|8]|9|) resembles very much 
the numerical simulation. The mobility above Walker 
breakdown is linear for the ID model curve, but shows a 
kink when the magnetization twist is taken into account. 
At the peak (hqH «38 mT), the velocity has increased 
from 4.7 m s -1 to about 10 m s" 1 , to be compared to 
the value of 9.2 m s _1 obtained in the simulation. At 
high fields, the calculated curve eventually coincides with 
the ID model one. The kink is in fact broader that the 
one observed in the simulation, which is not surprising 
considering the approximations needed to establish Eq. 
(8). A very similar result is obtained for the h=20 nm 
curve, showing that the velocity increase is due to the 
twist of the magnetization. 

A partial analytical integration of Eq. over t and z 
gives some insight as to which of the three terms H y , V 2 q 
and/or V 2- - absent from the ID model expression - gov- 
erns the shape of v(H). Taking into account the boundary 
conditions §f|o./i— ff lo./i— immediately shows that the 
V 2 ? and V 2 ip terms do not contribute at all. Using the 
odd- nature about z = h/2 of both H y and ip then leaves: 



A7 

< v >= 



1 + a- 



[ Ll M 1 f T T . „ . 27rf, . „ , , 7T f T . , f 

: }~ 2~ T J J o(2S sm— )sm20 o cft+ jj^ J sm ^ n J q 



siting / noH y (z) smtpdzdt + a hqH z > (10) 
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Jo(x) refers to the m=0 Bessel function of the first 
kind. In the absence of twist within the DW (<So=0), the 
ip term vanishes, and the ID expression is recovered. In 
the absence of time dependence of the twist (£(t)=So), 
the second term in Eg. (1101) vanishes by time integration 
of the T-periodic sin </>o term. This leaves an expression 
very similar to the ID equation, save for a decrease by 
Jo(2Sq) of the demagnetization term. 

When the twist is taken to be time dependent as de- 
scribed in Eqs. (3)- (5), the demagnetizing field (first term 
of Eq. (10)) created across the domain wall over one pe- 
riod is decreased compared to the pure Bloch DW case, 
as previously. The second term of Eq. (10) reflects the 
strength of the torque m x H y averaged over the thick- 
ness of the layer, and over a period. Because both the 
stray field and the z-dependence of the DW magnetiza- 
tion ip invert signs at mid height of the layer, the integra- 
tion over the layer thickness of this torque will always be 
non-zero. The sign of this contribution is moreover given 
by the phase of the twist amplitude S(t) with respect to 
the mean azimuthal angle of the DW <f>o(t). As shown in 
Fig. [3j when the twist is winding clockwise up the layer 
(S < 0), the magnetization lies opposite the direction of 
motion (00=0-180°, sin^o >0). When the twist is wind- 
ing counterclockwise (S > 0), the magnetization is facing 
the direction of motion (00=180-360°, sin</>o <0). As a 
result, this integration yields over a whole precession pe- 
riod a net positive contribution to the DW velocity. The 
amplitude of H y can easily be of the order of 1.3M at the 

surface, and so the m x H y torque will dominate the first 
term. Plotting only the applied field and demagnetizing 
field terms calculated with a twisted magnetization (la- 
beled A and B in the inset Fig. [5]) gives a curve very 
similar to the one computed from the ID model equation 
(6). It is thus clearly the stray field term (C) that yields 
a velocity increase. 

For an increased velocity to appear, it is therefore re- 
quired to combine both a twist of the DW magnetization 
across the layer, and a particular time dependence of this 
twist amplitude. As shown above, the twist and curva- 
ture are coupled variables. At fields giving a maximum 
flexion, a maximum twist will therefore also be obtained. 
The position of the resulting kink in the v(H) curve will 
then be indirectly related to the field dependence of the 
curvature tj(H), given by V 2 q. 

Finally, we come to the double kink for h—40 nm. At 
low field the kink is very similar to the /i=30 nm curve, 
but abruptly drops around mT. Since the veloc- 

ity increase is directly due to the continuous and symmet- 
rical twist of the DW across the layer's width, horizontal 
Bloch lines are likely to severely disrupt the scheme pre- 
sented above, all the more so as the HBL is tighter. As 
shown in Fig. HBLs precisely start appearing around 
fioH=30 mT for /i=40 nm. They are then possibly re- 
sponsible for splitting the peak into a sharp kink followed 
by a broad bump. 
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FIG. 5. [color on-line] Velocity computed analytically using 
the (q, tjj, Vq, V0) model for h~30 nm and H > Hw, and 
the corresponding simulated OOMMF curve, [inset] Close-up 
of the kink area evidencing the respective contributions to the 
velocity of the applied field (A), the demagnetizing field (B) 
and the stray field (C). See Eq. (9) and text for details. 



IV. POSITION OF VELOCITY PEAKS 
A. Simulations 

An important question remains concerning the posi- 
tion of the velocity peaks, as well as their thickness de- 
pendence. Because the shape of t]q(H) and Sq(H) recalls 
a resonance, the free oscillation of a DW was investi- 
gated. With H=0 in the numerical simulation, the DW 
is launched from one of the configurations shown in the 
top caption of Fig. QJ: either a "n=l"sine mode, or a 
"n=0" cosine mode, using once again the numbering of 
Slonczewski^I. The DW is then let to relax. In both 
cases, its movement is initially highly non-linear, and 
then rapidly falls into a damped oscillation of the n—1 
mode. Higher orders of the DW are also occasionally ob- 
served. The frequency spectrum of the DW elongation 
obtained by fast Fourier transform is centered around the 
free oscillation frequency Jfo 01 the DW. It is found to 
decrease when the thickness increases: fpo=8.20, 2.34, 
1.28 and 0.88 GHz for respectively h=10, 20, 30 and 
40 nm. 

Let us now compare these frequencies with the pre- 
cession frequency of the magnetization within the DW 
(solid line in Fig. [6j independent of layer thickness). As 
shown by the arrows in Fig. |6l when the applied field 
induces a precession frequency close to ]fo-, a velocity 
peak is obtained. The v(H) peak as well as the r)o(H) 
therefore likely correspond to the resonance of the n=l 
flexion mode of the DW. 

More exactly, the free oscillation resonance field falls 
slightly higher than the actual peak for h=20 nm (78 mT 
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instead of iioH res —70 mT) and 30 nm (48 mT instead of 
(ioH res =35 mT), and between the peak and the bump for 
/i=40 nm (35 mT). An explanation as to why the peaks 
are observed at lower frequency than expected from the 
free oscillation simulations can be put forward. A do- 
main wall in movement (under field) cannot be expected 
to have quite the same resonance frequency as a freely 
oscillating domain wall. In particular, applying an exter- 
nal field H z will disadvantage the higher order modes of 
the DW where magnetic charges due to its curvature will 
create an important magneto-static contribution. For 
H z = 0, the free oscillation frequency likely reflects the 
presence not only of the n—1 mode but also of higher 
order (and higher frequency) modes. These are indeed 
observed when plotting the DW depth-profile. These 
higher modes are progressively stifled upon increasing the 
field, and the frequency spectrum brought lower. This 
has been verified by doing Fig. [TJ;-type simulations for 
h=20 nm under an applied field of I, 3, 7.5 and 15 mT. 
As a result, the resonance under field appears at slightly 
lower frequency (and lower fields) than the free oscilla- 
tion frequency. An interesting conclusion to this is that 
a DW cannot truly be considered as a harmonic oscilla- 
tor as is often argued (Ref. and therein), but rather 
as a parametric oscillator whose resonance frequency will 
effectively depend on the forcing function (applied field 
in our case). 




100 200 300 

U H (mT) 

FIG. 6. Velocity versus field curves (closed symbols, shifted 
vertically for clarity), and precession frequency (solid line): 
the fields giving the peaks in the v(H) curves coincide with the 
ones yielding the simulated free oscillation frequencies (H=0). 
[inset] Comparison between the thickness dependence of the 
predicted resonance frequencies (/sl) given in Ref. [22| and 
the simulated ones (/fo)- 

Eq. (11) can serve as a general guide line, one must 
however keep in mind that the restoring force in Ref. 
[13] is somewhat different from the one in our geometry. 



B. Analytical model 

Slonczewski et ali^I established an approximate expres- 
sion for the flexural mode frequencies /sl of an assembly 
of domains arranged in stripes of infinite length (no ap- 
plied field), neglecting damping. In our 2D geometry, 
they can be expressed as a function of n-indexed wave- 
vectors k z (n) = mr/h: 



fsL = MoM + —k*^ MAk z + —kl (11) 

The inset of Fig. [6] compares the simulated resonance 
frequencies determined above (/fOi open squares), with 
the frequency of the first fsL mode computed from Eq. 
(11) with the numerical simulation's parameters. Since 
fsL was established at the first order in Afc z , it largely 
underestimates the resonance frequencies for very thin 
layers, but both plots eventually converge for h > 70 nm. 
Using Eq. (11) as a guideline shows that for h >> A and 
h > A, the position of the peaks in v(H) (given by Jfo) 
will vary roughly like 1/h with thickness and be very 
dependent on the A/h ratio. It will not depend on the 
uniaxial anisotropy. This is to be compared to the yA 
dependence of the oscillation frequency of a transverse 
DW simulated numerically in Permalloy^. For sample 
thicknesses much larger than the exchange length, the 
frequency will vary as l/Vh, A 1 / 4 and Ku 1 ^ ■ Although 



V. EXPERIMENTAL RESULTS 

This approach was then confronted to experimental 
data on perpendicularly magnetized samples. The di- 
lute magnetic semiconductors GaMnAs and GaMnAsP 
offer weak pinning which allows the observation of both 
the stationary and the precessional DW propagation 
regimes 1 ^. Moreover, magnetic parameters such as M, 
K u and A, can be adjusted by doping2i~— . Contrary 
to the virtual material used in the simulations, GaM- 
nAs also exhibits a complex in-plane anisotropy, which 
can induce anisotropic velocities leading to spectacular 
non-circular domain shapes^. However, the difference 
in velocity along the main axes being around 10 % only, 
the effect of the in-plane anisotropy is likely much weaker 
than the mechanisms described above. 

Fig. [7^ shows the v(H) curve at 20, 50 and 
80 K for a 50 nm thick Ga0.93Mn0.08As0.915P0.0s5 layer 
grown over GaAs. Fig. [7b,c,d show data for three 
Gao.93Mno.07As samples grown on a Inj,Gai_j,As pseudo- 
substrate: h=50 nm at temperatures 4-90 K (y=15%, 
previously published in Ref. [lfj), and h=20 nm (y=10%, 
T=40 K) and fc=40 nm (y=10%, T=60 K). After an- 
nealing, these last two layers have Curie temperatures of 
respectively 118 K and 136 K. The micromagnetic pa- 
rameters of both /i=50 nm layers were determined using 
Kerr microscopy as detailed in Ref. [2l[ (where they are 
referred to as samples B and D). The main differences be- 
tween these two samples are a higher uniaxial anisotropy 
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FIG. 7. [color on-line] Experimental velocity versus field data 
(curves shifted for clarity where needed). The polygons indi- 
cate the velocity anomaly regions. The arrows point to the 
H=0 resonance fields determined numerically and the cir- 
cles to the lower parametric excitation field (see text), (a) 
Ga0.93Mn0.08As0.0915P0.0s5, h=50 nm. (b) Gao.93Mno.07 As, 
h=50 nm (reproduced from Ref. [lfl])- (c) Gao.93Mno.07As 
h=20 nm. (T=40 K), (d) and h=40 nm (T=60 K). 



in GaMnAs, and a higher exchange constant and satu- 
ration magnetization in GaMnAsP. In all these measure- 
ments, the DW velocity was determined using magneto- 
optical imaging and a magnetic field pulse technique fully 
described in Ref. [Toj [ . 

All the v(H) curves evidence low-field and high-field 
anomalies to the ID model curve broadly indicated by 
polygons in Fig. [7] In the GaMnAs data, the /i=50 nm 
curves show a weakly temperature-dependent kink at 
around 100-110 mT, and a discreet one at lower field 
(55 mT at T=4 K). The latter rapidly merges with the 
Walker peak ([IqHw ~5 mT) when temperature is in- 
creased. At high temperature, a negative mobility region 
is eventually obtained. The /i=40 nm curve is similar: a 
broad kink around Walker break-down (imqHw ~5 mT), 
followed by a "split" bump centered around 100 mT. 
The h=2Q nm curve exhibits distinctly Walker break- 
down (/zo-HV = 14 mT), a narrow negative mobility region 
before a peak at 23 mT. It is followed by a broad split 
bump centered around 105 mT. The high field mobility 
is similar for h=A0 nm and h=50 nm, and double that of 
the h=20 nm layer. To summarize, all GaMnAs samples 
exhibit a first velocity anomaly varying with temperature 
and thickness, and a second one independent of those pa- 
rameters at «105 mT. The simulations (h=20 nm and 
40 nm of Fig. Hb,d) reproduce qualitatively the low-field 
anomaly of the corresponding GaMnAs samples, but not 
the high-field 105 mT one. The GaMnAsP data (Fig. 
[TJi) exhibits similarly two anomalous velocity features: a 
zero or negative mobility regime at low field and a very 
broad bump appearing at higher field: 50 mT at T=20, 



50 K and 34 mT at T=80 K. 

The resonance frequencies were then determined for 
each /i=50 nm sample by micromagnetic simulations 
(fioH—0) as described above, using the experimental val- 
ues of M, K u and A^. No sufficiently accurate value of 
A was available for the ft.=20, 40 nm samples. The res- 
onance fields found by comparing fpo to the precession 
frequency are indicated by arrows directly on the v(H) 
curves of Fig. [JJi,b: 48, 40 and 23 mT for the GaMnAsP 
sample (T=20, 50 and 80 K), and 29 mT for the GaMnAs 
sample at T=80 K. Since the simulations showed that the 
low-field velocity anomalies appeared slightly lower than 
the resonance fields, the arrows very likely correspond to 
the flexion/twist-induced kinks, in good agreement with 
the description of the shifted peak position developed in 
Section IV. This would explain qualitatively the fact that 
the negative mobility regime appears upon increasing the 
temperature, as the flexion-induced kinks down-shift to- 
wards the Walker field. Taking into account the probable 
exchange constant inhomogeneity across the layer thick- 
ness would likely give better quantitative description of 
the kinks. 

The «105 mT kinks in GaMnAs and «60 mT broad 
bumps in GaMnAsP can therefore not be explained by 
the twist /curvature model, as they shift only weakly 
with thickness or temperature. A tentative explana- 
tion for the origin of these kinks is the parametric ex- 
citation put forward by Randoshkin 35 whereby the bulk 
spin waves of the domains couple to the DW magne- 
tization, thus providing an efficient energy dissipation 
channel leading to a velocity increase. This phenomenon 
is expected for frequencies equal to twice the preces- 
sion frequency: fbuik—^fprec- This corresponds to res- 
onance fields around H u /3 or 2H u /3 depending whether 
the bulk spin waves are emitted in the domain ori- 
ented anti-parallel or parallel to the applied field re- 
spectively, and using the usual expression of bulk spin 
waves frequencies 3 ^. When anisotropy coefficients were 
available, the lower bound expression could be evalu- 
ated (circles on the curves of Fig. [7^,,b,c). The higher 
bound term falls out of the explored field range. For the 
GaMnAs samples, it yields about 110-140 mT (T=12- 
100 K) for h=50 nm, wl40 mT for h=20 nm, and for 
GaMnAsP 61-78 mT (T=20-80 K) . These anharmonic 
coupling fields follow reasonably the temperature depen- 
dence of the high-field bumps observed for both GaM- 
nAs and GaMnAsP, but do not explain the broadness of 
the bumps. The mechanism for this coupling has more- 
over never been fully investigated. In the simulations, it 
would appear at around H u /3=179 mT. It has not been 
observed, possibly due to the strong damping in the sim- 
ulations. 



VI. CONCLUSION 

To summarize, we have investigated numerically field- 
driven DW propagation in ferromagnetic layers of thick- 
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ness close to the exchange length. These simulations re- 
veal anomalous velocity peaks appearing in the preces- 
sional regime. They are correlated with the maximum 
flexion of the DW and the maximum twist of the magne- 
tization vector inside the wall, both being intimately cou- 
pled. In order to elucidate their respective roles, we com- 
plement the historical (q, ip) DW propagation model with 
the phenomenologically determined VV> (twist). The ve- 
locity peaks are then well reproduced and found to result 
from the torque effect of the stray field on the twisted 
magnetization. When the DW magnetization is suffi- 
ciently twisted across the layer thickness, the stray field 
emanating from the up- and down-domains surrounding 
it induces an efficient torque on the total DW magnetiza- 
tion. Because the amplitude of this twist varies over one 
precession period, the torque averages out to a strictly 
positive value. Anomalies in v(H) curves are therefore 
obtained for the particular range of thickness over ex- 
change length ratio allowing a strong twist of the DW 
magnetization, but forbidding the appearance of horizon- 
tal Bloch lines, which tend to stifle the velocity enhance- 
ment. The resonance field of the twist/curvature is such 
that the corresponding precession frequency is then close 
to the DW fiexural resonance frequency, and strongly 
depends on the layer thickness. The experimental v(H) 
curves obtained for ferromagnetic GaMnAs and GaM- 



nAsP layers show velocity enhancement related to this 
DW magnetization twist and stray field-induced torque 
just above Walker break-down. At high field, they exhibit 
broad velocity bumps of different origin, which may be 
related to non-linear excitation of bulk spin-waves not 
evidenced by our 2D simulations. A full 3D simulation 
of the velocity would surely yield additional information. 

The description of the velocity versus field curve in 
terms of the (q,V',Vq,V'!/ , ) variables is an alternative 
and promising way to understand DW propagation ex- 
periments in perpendicular-to-plane magnetized samples. 
Keeping in mind that the case of Bloch DWs is sim- 
pler in essence, the general ideas of this model also 
seem applicable to in-plane magnetized materials. Ve- 
locity versus field curves taken in Permalloy nanostrips 
have for instance very clearly evidenced kinksi^ in the 
precessional regime at fields decreasing with nanostrip 
width, the relevant dimension for the flexion when m is 
in the plane of the sample. This would however call for 
more thorough investigation, as well as the relevance of 
the exposed model in current-induced DW propagation 
experiments^. 

We acknowledge L. Largeau for the TEM determina- 
tion of the GaMnAsP layer thickness. This work was in 
parts supported by Region He de France under contract 
IF07-800/R with C'Nano IdF. 
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